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ABSTRACT 

Using the Gemini South 8m telescope, we obtained high resolution 
11.7 and 18.3 /im mid-IR images of SN 1987A on day 6526 since the ex- 
plosion. All the emission arises from the equatorial ring. Nearly contem- 
poraneous spectra obtained at 5-38 /im with the Spitzer Space Telescope 
show that this is thermal emission from silicate dust that condensed out 
in the red giant wind of the progenitor star. The dust temperature is 
166^^2 K; an d the emitting dust mass is (2.6j£2) x 1(T 6 M . Compar- 
ison of the Gemini 11.7 /im image with Chandra X-ray images, Hubble 
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UV- optical images, and ATCA radio synchrotron images shows gener- 
ally good correlation across all wavelengths. If the dust resides in the 
diffuse X-ray emitting gas then it is collisionally heated. The IR emission 
can then be used to derive the plasma temperature and density, which 
were found to be in good agreement with those inferred from the X-rays. 
Alternatively, the dust could reside in the dense UV-optical knots and be 
heated by the radiative shocks that are propagating through the knots. In 
either case the dust-to-gas mass ratio in the CSM around the supernova 
is significantly lower than that in the general interstellar medium of the 
LMC, suggesting either a low condensation efficiency in the wind of the 
progenitor star, or the efficient destruction of the dust by the SN blast 
wave. Overall, we are witnessing the interaction of the SN blast wave 
with its surrounding medium, creating an environment that is rapidly 
evolving at all wavelengths. 

Subject headings: Stars: Supernovae: Individual: SN 1987A — Infrared: 
ISM: Dust, Supernova Remnants 



1. INTRODUCTION 

Since its explosion, SN 1987A has evolved from a supernova (SN) dominated 
by the emission from the radioactive decay of 56 Co, 57 Co and 44 Ti in the ejecta to a 
supernova remnant whose emission is dominated by the interaction of the supernova 
blast wave with its surrounding medium. The medium surrounding the SN is domi- 
nated by the well-known "circumstellar envelope" (CSE), which consists of an inner 
equatorial ring (ER) flanked by two outer rings (Burrows et al. 1995), possibly part 
of an hour-glass structure. 



1 Based on observations obtained at the Gemini Observatory, which is operated by the Associa- 
tion of Universities for Research in Astronomy (AURA), Inc. under cooperative agreement with the 
NSF on behalf of the Gemini partnership: the National Science Foundation (United States), the 
Particle Physics and Astronomy research Council (United Kingdom), the National Research Coun- 
cil (Canada), CONICYT (Chile), the Australian Research Council (Australia), CNPq (Brazil), and 
CONICET (Argentina). 
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The collision between the ejecta of SN 1987A and the ER predicted to occur 
sometime in the interval 1995-2007 (Gaensler et al. 1997; Borkowski et al. 1997) is 
now underway. At UV-optical (UVO) wavelengths, "hot spots" have appeared inside 
the ER (Pun et al. 1997), and their brightness varies on time scales of a few months 
(Lawrence et al. 2000). New hot spots continue to appear as the whole inner rim 
of the ER lights up. The visible-light HST image obtained in 2004 Feb 20 reveals 
a necklace of such hot spots which nearly fill a lighted ring. Ongoing monitoring at 
X-rays with the Chandra and at radio frequencies, shows that the evolution of the 
emission from the ER follows a similar pattern at all wavelengths. 

There exist very few mid-infrared (IR) observations of supernovae in general. 
Therefore SN 1987A, the closest known supernova in 400 years, gives us an oppor- 
tunity to explore the mid-IR properties of supernovae and the dust in their ejecta 
and surrounding medium with the help of the newest generation of large-aperture 
telescopes and sensitive mid-IR instrumentation like the T-ReCS on the Gemini, in 
combination with IR data obtained from Spitzer Space Teles cope(S ST) (Werner et 
al. 2004). The T-ReCS observations of the mid-IR emission from SNR 1987A are 
part of our continuous monitoring of the SN and its surrounding medium. The first 
detection and analysis of mid-IR emission at the position of the supernova has been 
reported in (Bouchet et al. 2004) (hereafter Paper I). 

The origin of the mid-IR emission could be line emission from atomic species, 
synchrotron or free-free continua, or thermal emission from dust which is proba- 
bly the dominant source of emission. In general, there are several scenarios for the 
origin and the heating mechanism of the dust giving rise to the late time mid-IR 
emission in Type II supernovae (Graham & Meikle 1986; Gerardy et al. 2000). 
Thermal mid-IR emission could be: (1) the emission from SN-condensed dust that is 
collisionally heated by reverse shocks traveling through the SN ejecta; (2) the emis- 
sion from circumstellar /interstellar dust heated by the interaction of the expanding 
SN blast wave with the ambient medium; and (3) the delayed emission (echo) from 
circumstellar dust radiatively heated by the early UVO supernova light curve. 

Our imaging observations provide strong constraints on these possible scenarios 
for the mid-IR emission. They show that the bulk of the mid-IR emission is not 
concentrated on the center of the explosion, but arises from a ring around the SN. 
The morphology of the mid-IR emission can therefore be used to eliminate several 
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scenarios for its origin. Combined with Spitzer spectroscopy, our observations can be 
used to determine the dust composition and temperature distribution in the ring. The 
Chandra X-ray, and HST UVO data provide important constraints on the physical 
conditions of the medium overtaken by the supernova blast wave. These constraints 
can be used to determine the physical association and the heating mechanism of the 
dust giving rise to the mid-IR emission. 

The paper is organized as follows: We first describe in §2 the imaging observa- 
tions of the SN obtained by T-ReCS, and compare its IR morphology to that at radio, 
X-ray, and UVO wavelengths. Lower resolution Spitzer mid-IR imaging observations 
detected SNR 1987A (including the ring) as an unresolved point source, and are 
used in conjunction with spectroscopic data to determine the possible contribution 
of lines, and the composition of the dust giving rise to the continuum emission. In 
§3 we describe the procedure used for the analysis of the spectroscopic and imaging 
data, determining the dust composition, and presenting maps of dust temperature, 
IR opacities, and dust column densities. The IR image of the circumstellar medium 
around SN 1987A has a morphology similar to the Chandra X-ray and HST UVO im- 
ages. The limited mid-IR resolution does not allow us to unambiguously determine 
whether the dust resides in the X-ray emitting gas, or in the UV-optical line emitting 
knots in the ER. We therefore resort in §4 to an analysis of possible dust heating 
mechanisms and to calculating the inferred dust masses and dust-to-gas mass ratios 
for several possible scenarios. In §5 we discuss the evolution of the supernova and 
its environment as manifested from the observed light curves at various wavelengths. 
The results of our paper are summarized in §6. 



2. OBSERVATIONS 

2.1. Mid-Infrared Gemini Observations 

The T-ReCS mid-IR imager /spectrometer at the Gemini 8m telescope offers a 
combined telescope and instrument with diffraction limited imaging (~ 0.3" resolu- 
tion) and superbly low thermal emissivity. On 2003 Oct 4 (day 6067), we imaged 
SNR 1987A with T-ReCS as part of the instrument's System Verification program 
and we reported on the detection of 10 and 20 /zm emission from the ER, and on a 
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10 /im emission from the supernova's ejecta (Paper I). Subsequent observations were 
carried out in January 6, 2005 (day 6526) in the narrow Si5 filter (X e ff = 11-66 /im; 
> 50% transmission at A = 11.09-12.22 /im), and in February 1, 2005 (day 6552) 
in the Qa filter (X e ff = 18.30 /im; > 50% transmission at A = 17.57-19.08 /an). 
Results are presented in Figures 1-a and 1-b. These images show several luminous 
"hot" spots distributed over the ring. The calibrated flux density integrated within 
an aperture of 1.3 arcsec radius is F u (11.7 /im) = 18.4 ±1.2 mJy in the Si5 filter, 
and ^(18.3 /im) = 53.4 ±9 mJy in the Qa filter. No color correction was applied 
and this would most likely increase the flux density. The standard star used for the 
calibration of the 11.7 //m measurement was HD 29291, whose flux density was taken 
to be 6.78 Jy at 11.7 /im. We used a CMa with a flux density of 44.3 Jy at 18.30 
/im for the flux calibration of the 18.3 /im observation. 

The black body colour temperature corresponding to the measured fluxes at 
these two wavelengths is T = 185 K and the luminosity is L B b = 3.74 x 10 36 erg s -1 . 
We use the Mathis (1990) extinction law with t 18 .3 = Tn.y/1.35 and A 18 _ 3 /Aj = 0.083 
and Au_ 7 /Aj = 0.098 to compute the black body temperature and the optical depth 
for each individual pixel, resulting in the maps shown in Fig 1-c and 1-d. Note that 
in order to have the iterative algorithm converge we must assume reasonable values 
as a starting point. We stress that Fig 1-c and 1-d show the color temperature 
and optical depth maps, which are slightly different from the maps related to the 
physical dust as calculated in Section 3. They are shown only for illustrating the 
results obtained from our data fitted to the simplest modeling (eg. black body). 

In Figure 2 we compare our 11.7 /im new data with the one obtained in the 
broad N band (10 /im) on October 2003 at day 6067 (Paper I). This figure shows 
a clear brightening in the South- West region of the ER, superimposed on a general 
brightening all over the ring. Given that the two images are taken with significantly 
different filters we investigated whether the difference could be due to the different 
spectral coverage: the spectra obtained with Spitzer which are discussed in next 
section do not show any feature that could be a source of the difference. Thus we 
conclude that we are observing a true brightness change. 

Figure 3(a,b) displays our images in both filters with the contours of the 0.3-8 
keV X-ray image from the Chandra Observatory obtained nearly simultaneously at 
day 6529 (January 9-13, 2005) (Park et al. 2005b). The correspondences between 
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our 11.7 jum and that obtained by the Chandra is very good, but less so at 18.3 /im 
because of the lower signal-to-noise ratio in that image. A more detailed discussion 
of the relation between the IR and X-ray images of the ER will be presented in §3 
below. 

Figure 4(a, b) shows the contours of the image obtained in the 12 mm band (16- 
26 GHz) at day 6003 (July 31, 2003) with the Australian Telescope Compact Array 
(ATCA) at the Australian National Telescope Facility (ATNF) (Manchester et al. 
2005) superimposed on our 11.7 and 18.3 /im T-ReCS images, respectively. The 
correspondences between our 11.7 /im image with the synchrotron radio emission is 
not as remarkable as it is with the Chandra image. Nevertheless, the 18.3 /im bright 
spot in the East side looks better correlated with the radio lobe than X-ray spots do. 
Furthermore, an image obtained in May 5, 2004 (day 6298) at the same frequencies 
is posted on the ATNF web page and is reproduced here in Figure 4(c,d). This 
later radio image shows better correlation with our 11.7 /im image than the radio 
image obtained on day 6003, probably because it was taken closer to the epoch of 
the mid-IR observations. This demonstrates the importance of evolutionary effects 
on the morphology of the emission at all wavelengths. 

There is good overall agreement in shape and size between our IR image and 
images obtained from the X-ray to the radio. The mean radii and approximate sur- 
face brightness distribution (brighter on the east side) of the ring are similar at all 
wavelengths, demonstrating that the dust is co-extensive with the gas components. 
The origin of that brightness asymmetry may be related to an asymmetric distribu- 
tion of the ejecta or the CSM (Park et al. 2002, 2004), and/or to a time- dependence 
effect caused by the tilt of the ER as argued by Panagia et al. (1991). 

The most likely source for mid-IR radiation is thermal emission from warm dust 
(see discussion below). The X-ray radiation is thermal emission from very hot gas 
(Park et al. (2005a), and previous references therein), whereas the optical emission 
arises from the dense knots in the ER that are overrun by slower shocks. The 
radio emission is likely to be synchrotron radiation from shock-accelerated electrons 
spiralling in the remnant's magnetic field as stated by Dunne et al. (2003) and 
Manchester et al. (2005). Park et al. (2002, 2004) argue that, until 2000 December, 
hard (E > 1.2 keV) X-ray and radio emissions were produced by fast shocks in the 
CS HII region while the optical and soft (E < 1.2 keV) X-ray emissions came from 
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slower shocks in the denser ER. Park et al. (2003) note that as of 2002 Dec. 31 (day 
5791) correlations between the X-ray and the optical/radio images are more complex 
than the above simple picture, which is expected as the blast wave is reaching the 
main body of the inner ring. 

2.2. Spitzer Observations 

2.2.1. Imaging Data 

Imaging of SNR 1987A was carried out with Spitzer Space Telescope's MIPS 
instrument at 24 /iin (Rieke et al. 2004) and IRAC instrument at 3.6 - 8 /im 
(Fazio et al. 2004) (AORIDs = 5031424 and 5030912). Almost one year later 
SNR 1987A was again imaged with IRAC, serendipitously near the edge of the field 
of observations targeting other sources (AORIDs = 11191808 and 11526400). All 
these data were obtained from the Spitzer data archive. The corresponding images 
are shown in Figure 17(a-f), which shows also a near-IR image from the Hubble 
Space Telescope for comparison. At 24 and 8 /im SNR 1987A was detected as an 
unresolved point source amid a field of complex cirrus emission. At 5.8 /im, the 
source appears very slightly distorted, and at 4.5 and 3.6 /im the SN appears to be 
swamped by the emission of companion stars 2 and 3. The flux densities at 24, 8, 
and 5.8 /im were measured using SExtractor (Bertin & Arnouts 1996) to perform 
aperture photometry on the post-BCD images. Aperture radii used were 4, 5, and 
6 pixels (4.8", 6.0", and 14.7") at 5.8, 8.0, and 24 /im respectively. Approximate 
aperture corrections of 1.10, 1.07, and 1.14 were applied using information from the 
IRAC and MIPS Data Handbooks. The calibrated flux densities are given in Tab. 1. 

2.2.2. Spectroscopic Data 

Observations of SNR 1987A were also performed with Spitzer^ Infrared Spec- 
trograph (IRS) in its short (wavelength) - low (resolution), short-high, and long-high 
modes (Houck et al. 2004). These data were also obtained from the Spitzer data 
archive (AORID = 5031168). Figure 18 shows the slit positions for the different sets 
of observations. For each spectral mode we extracted SN 1987A spectra from the 2-D 



- 8- 



coadded post-BCD images using SPICE (http://ssc.spitzer.caltech.edu/postbcd/spice.html). 
For the short-low data, the source is placed at 4 positions along the slit to generate 
spectra at two positions for two different spectral orders. For each slit position, a 
2-D background image was generated from the median value of the images for the 
other three slit positions. These backgrounds were subtracted prior to extracting 
spectra with SPICE. For the high resolution observations, only two slit positions are 
observed within a much narrower slit. So for these observations, SPICE was used 
to extract the SN spectrum from the columns occupied by the unresolved source, 
and a background spectrum from columns toward the opposite end of the slit. The 
background spectra were then subtracted from the source spectra. For all modes 
the spectra from the two slit positions (per spectral order) were combined using a 
weighted (l/cr 2 ) average, and (generally noisy) data where spectral orders overlapped 
were discarded. For the short-hi data, an empirically determined scaling factor of 
1.38 was applied before averaging to bring results from both slit positions into agree- 
ment. An additional scale factors of 1.46 was applied to the short-high spectrum to 
normalize it to the short-low data, and a subsequent factor of 1.25 was applied to 
long-high data to normalize to the short-high data. These latter scaling factors are to 
be expected if the subtracted background was partially contaminatead by the source, 
or in the case that the smaller high resolution slits were slightly misaligned. Finally 
for the sake of comparison with the broad band measurements and dust models, the 
high resolution spectra were median-binned using intervals of 12 and 1 1 wavelength 
samples for the short and long wavelengths respectively. 

Figure 5 shows the overall calibrated spectrum and Figure 6 displays individual 
profiles and identifications of the main emission lines detected with IRS from the 
Spitzer Space Telescope: it shows that the T-ReCS observations are dominated by 
the dust continuum emission, and not the lines. The [Ne II] 12.81 /zm and [Ne III] 
15.56 /iin lines are however clearly seen. A weak [Si II] 34.81 /xm line remains after 
background subtraction, while [S III] 33.48 /zm disappears entirely. Two strong lines 
are seen near 26 /xm: the redder line is [Fe II] at 25.99 /im and the bluer line is [O 
IV] at 25.89 /im. Both lines could be arising from [Fe II] (or both [O IV]) if from fast 
moving ejecta on near and far sides of the explosion, but the lack of splitting of the 
other lines makes this seem unlikely. A weak third component seems present here as 
well. If fitted as an unresolved line this may be [F IV] at 25.83 fxm. However, it may 
also be fitted as the wing of a very broad line (FWHM ~ 2800 km s _1 ) underlying the 
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narrower [O IV] and [Fe II] lines. Such high velocity would indicate an association 
with the SN ejecta rather than the ER for this line. Line fluxes, centers, and widths 
have been calculated by gaussian fits using SMART (Higdon et al. 2004) and results 
are given in Tab. 2. 

2.3. Near-IR CTIO Observations 

Near-IR J(1.25 pm), H(1.65 fjm), and K(2.2 /im) imaging observations of SNR 1987A 
were obtained on 2005, January 3-5 with ISPI attached to the Blanco 4-m telescope 
at the Cerro Tololo Interamericain Observatory. Results are displayed in Figure 7. 
We clearly see the ER at J and K, while most of the emission detected in the H band 
arises from the supernova itself. These images have been deconvolved by using the 
"Multiscale Maximum Entropy Method" algorithm developed by Pantin & Starck 
(1996). The flux densities measured in the CTIO data are (Fj = 0.33 ± 0.25 mJy; 
F H = 0.11 ± 0.05 mJy; F K = 0.51 ± 0.2 mJy). These fluxes refer to the total of 
the ER + ejecta in each band (although each band seems dominated by one or the 
other). Images obtained with HST on November 11, 2005 through the filters F110W, 
F160W, and F205W, are also shown in this figure, and show that compared to the 
ring, the ejecta are relatively brighter in H band (F160W) than in other bands. This 
clearly demonstrates the validity and the power of the deconvolution algorithm ap- 
plied to our CTIO data. Note, however, that the HST bands are not exactly J, H, 
K: one significant difference is that the Paschen a line is missed between the ground 
based H and K bands, but would be seen in HST's F205W band. 

The Spitzer mid-IR spectrum indicates that dust in SN 1987A and the ER is 
too cool (see Section 3) to emit significantly in these near-IR bands. We believe that 
line emission dominates these broadband near-IR observations, and that differences 
in the composition, temperature, and density of the ejecta vs. CSM lead to the 
different relative brightnesses of these structures. The young supernova remnants 
Cas A and Kepler may be the closest analogs we have for interpreting SN 1987A's 
near-IR emission. Near-IR spectra and imaging of these SNRs by Gerardy & Fesen 
(2001) and Rho et al. (2003) reveal similar variation in near-IR colors. Circumstellar 
material (e.g. Cas A's quasi-stationary flocculi) is bright in J due to He I 1.083 /im, 
and emits strong [Fe II] lines in both J and H bands. Weaker hydrogen lines can be 
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detected in all bands. The ejecta in several of Cas A's fast-moving knots emit strong 
[S II] lines at 1.03 /im in J, and weaker lines in [Fe II] and Si in various ionization 
states. Thus, the near-IR emission of the circumstellar ER may be dominated by 
H and He I lines in J and H with synchrotron emission possibly augmenting the K 
band (Gerardy & Fesen 2001; Rho et al. 2003). The synchrotron processes deserve 
consideration in SN 1987A particularly in view of the fact shown in Figure 4 that 
the radio morphology of SN 1987A reproduces the ring shape now. Also, the H2 
molecule has many emission lines in the K-band and longer wavelengths. Spectra 
are clearly needed for testing the plausibility of the contribution from this molecule 
or otherwise. The first overtone band head of molecular CO lies at 2.29/im and may 
be also contributing to K-band emission in the ring. The relatively strong H band 
emission of SN 1987A's ejecta could reflect [Fe II] emission, because a large fraction 
of the ejecta should be Fe from the core of the explosion, whereas Cas A's FMKs 
generally consist of lighter metals originating in the outer layers of the progenitor. 



3. DATA ANALYSIS 

3.1. The Dust Properties 

Figure 5 clearly demonstrates that the mid-IR emission is dominated by thermal 
emission from dust. The specific luminosity of a single dust particle of radius a at 
temperature T d , at wavelength A is given by: 

4(A) = Arra 2 nB u (X, T d )Q{\) (1) 
= 4m d «(A, a)nB v {\, T d ) 

where B U (X, T d ) is the Planck function, Q(X) the dust emissivity at wavelength A, 
and k(X, a) = 3Q(X)/4pa is the dust mass absorption coefficient, where p is the mass 
density of the dust particle. In the Rayleigh limit, when a < A, k is independent of 
particle radius. Figure 8 illustrates the k(X) curves for amorphous carbon (Rouleau 
& Martin 1991), graphite and silicate grains (Laor & Draine 1993) over the 5 to 
30 /im wavelength region. The figure shows the distinct optical properties of these 
dust particles over the mid-IR wavelength regime which can greatly facilitate the 
identification of the emitting material with even limited broad band filters. 
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For an optically thin point source, the flux density, F V (X), at wavelength A is 
given by: 

FJX) = 4 Md m ^(a. t A (2) 

where M d is the dust mass, and D is the distance to the supernova, taken to be 
D = 51.4 kpc (Panagia 1999). 

For an extended optically thin source with an angular size Q, the surface bright- 
ness, i,/(A), is given by: 

J„(A) = n T d (\)B v (\, T d ) (3) 
where T<j(A), the dust optical depth, is given by: 

M d k{\) 

r d {\) = (4) 



3.2. Spectral Analysis 

We fitted the integrated T-ReCS flux densities with a population of dust particles 
consisting of a single population of bare graphite, silicate or amorphous carbon grains, 
using equation (2). Optical properties of silicate and graphite grains were taken 
from Laor & Draine (1993), and those for amorphous carbon (BE) were taken 
from Rouleau & Martin (1991). The results are given in Figure 9, which shows 
that the T-ReCS observations alone cannot discriminate between the different dust 
compositions. Comparison with the Spitzer observations clearly demonstrates that 
the IRS observations can be well fit with a silicate dust composition (Figure 10), 
ruling out graphite or carbon dust as major dust constituents in the CSM. 

Figure 10 shows the fit of the silicate dust spectrum to the Spitzer IRS spectrum. 
The silicate temperature is T d = 180^^ K, and the dust mass is M d = (1.11q|) x 
1O~ 6 M . The global parameters resulting from our model for silicates are given in 
Table 3. Figure 11 shows the residuals of the fit, obtained by subtracting the silicate 
model fit from the data. The residuals are small, and their spectrum is too sharply 
peaked at short wavelengths to be fitted with any blackbody. The residuals are also 
broader than typical atomic line widths, suggesting that they may be due to solid 
state features in the dust, reflecting differences in the crystalline structure of the 
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silicates in the CSM from the average interstellar silicate dust used in the model. 
This figure shows also that the residuals cannot be fit by emission from amorphous 
carbon dust (shown by the green line in this figure). 

3.3. Image Analysis 

Our spectral analysis showed that the mid-IR spectrum of the ER is dominated 
by silicate emission. We therefore used the 11.7 and 18.3 fxm images of the remnant 
to construct temperature, and dust opacity maps of the circumstellar ring using Eq. 
3 to derive the dust temperature, and Eq. 4 to derive the dust optical depth. In 
calculating these quantities, we applied a background threshold of 0.03 mJy/pix at 
11.7 /im and of 0.08 mJy/pix at 18.3 /zm. The mass of the ring was calculated from 
the optical depth, over a surface area of 269 pixels, corresponding to the number of 
pixels that had a flux exceeding the respective thresholds at each wavelength. The 
average dust temperature in these maps is T d = I66TJ12 K. The average 11.7 /im 
optical depth per pixel is r d = (5. ^L^Jf) x 10~ 6 , and the total dust luminosity is 
L d = (2.3±S:t) x 10 36 erg s _1 , giving a dust mass M d = (2.6±?;2) x 10~ 6 M , in 
good agreement with the total mass obtained from the spectral analysis of the ER. 
Figure 12 shows the maps of the silicate dust temperature (a) and optical depth (b) 
in the ER. 



4. THE ORIGIN OF THE MID-INFRARED EMISSION 

The data presented in this paper show unambiguously that the emission is ther- 
mal emission from dust. At issue are the location and heating mechanism of the 
dust. The forward expanding non-radiative blast wave is currently interacting with 
the circumstellar material and the knots in the ER. The interaction of the blast wave 
with the knots transmits lower velocity radiative shocks into these dense regions, pro- 
ducing soft X-rays and the "hot spots" seen in the HST images. The interaction of 
the blast wave with the dense knots also generates reflected shocks that propagate 
back into the medium that was previously shocked by the expanding SN blast wave. 
The complex morphology and density structure of the ER gives rise to a multitude 
of shocks characterized by different velocity, temperatures, and post shock densities. 
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The mid-IR images cannot determine the location of the radiating dust, whether 
it resides in the X-ray emitting gas or in the denser UVO emitting knots. Therefore, 
we can not, a priori, assume a particular dust heating mechanism: collisional heating 
in the shocked gas, or radiative heating in the radiative shocks. 



4.1. Dust Heating Mechanism 



The relative importance of the two dust heating mechanisms is given by the 
ratio (Arendt, Dwek, & Moseley 1999): 

^ _ H rad _ n e n H A(T e ) P abs 
H coil n d n e A d (T e ) 



where H rad and H co u are, respectively, the radiative and collisional heating rates of 
the dust, n e , T e , and A(T e ) are, respectively, the electron density, temperature, and 
the atomic cooling function (erg cm 3 s _1 ) of the gas, P a bs is the fraction of the cooling 
radiation that is absorbed by the dust, and A d (T e ) is the cooling function of the gas 
via electronic collisions with the dust and given by: 

A d (T e ) = 2v e ira 2 kT e (h) (6) 

where n d the number density of dust grains, a their average radius, v e = (SkTz/ixme) 1 / 2 
is the mean thermal speed of the electrons, and (h) < 1 is the collisional heating 
efficiency of the dust which measures the fractional energy of the electrons that is 
deposited in the dust (Dwek et al. 1987a). 

In the following we will examine possible locations of the dust giving rise to the 
IR emission. For each possible site, the X-ray emitting gas or the optical knot, we 
determine the dominant cooling mechanism, the temperature of the dust, and the 
inferred dust mass. A site is viable if it can maintain a dust temperature between ~ 
150 and 200 K, the range of values reflecting the range of dust temperatures derived 
from the T-ReCS and Spitzer observations, and if the derived dust mass does not 
violate any reasonable abundances constraints. 
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4.2. Dust in the X-ray Emitting Gas 

The morphological similarities between the 11.7 //m emission, dust temperature, 
and optical depth maps on one hand and the X-ray maps of the supernova on the 
other hand suggests that the dust giving rise to the IR emission may be well mixed 
with the X-ray emitting gas. 

For an optically thin plasma P a b s ~ = ndira 2 (Q)£, where (Q) is the radiative 
absorption efficiency of the dust averaged over grain sizes and the radiation spectrum, 
and i is a typical dimension of the emitting region. Inserting this value in eq. (5) 
we get: 

n e A(T e )i <Q> m 
' " 2v e kT e (h) { ) 

7im e n e A(T e ) I (Q) 

X — 71 , ~ — X 



32 (A;T e )3/2 ( h ) 

For average conditions in the X-ray plasma (see below), characterized by plasma 
temperatures and densities of T e ps 10 7 K, and n e 300 — 10 3 cm -3 , we get A(T e ) rs 
4 x 10" 23 erg cm 3 s -1 , and (/i) > 0.1 for dust particles with radii larger than 0.05 jum. 
Adopting a value of (Q) ~ 1, gives an upper limit of 

K < 2.1 x 10~ 20 £ (8) 

The size of the X-ray emitting plasma is less than the radius of the ER which is £ < 
0.7 lyr = 7 x 10 17 cm, giving 1Z < 0.14, that is, the dust heating in the X-ray gas is 
dominated by electronic collisions. 



4-2.1. The Dust Temperature 

Figure 13 depicts contours of the temperature of collisionally heated silicate 
grains of radius a = 0.10 fj,m as a function of plasma density and temperature. 
To calculate the energy deposited in the dust we used the electron ranges of Iskef, 
Cunningham, & Watt (1983) for energies between 20 eV and 10 keV, and those 
of Tabata, Ito, & Okabe (1972) for higher incident electron energies. We also 
assumed that all incident electrons penetrate the dust (no reflection). Above gas 
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temperatures of about 3 x 10 6 K, the temperature of collisionally heated dust is 
independent of grain radius, and very well represented by a single dust temperature, 
since both the radiative cooling and the collisional heating rates are proportional to 
the mass of the radiating dust particle. Furthermore, at these plasma temperatures 
the dust temperature is essentially independent of gas temperature as well, and 
therefore an excellent diagnostic of plasma densities. The figure shows that for 
plasma temperatures above ~ 3 x 10 6 K, dust temperatures between 150 and 200 K 
require electron densities of about 300 to 1400 cm~ 3 . 

X-ray observations show the presence of two X-ray emission components (Park 
et al. 2006): one associated with the "slow shock" with an electron temperature of 
~ 0.23 keV and a density of n e ~ 6000 cm" 3 , and the second, associated with the 
"fast shock" with electron temperatures and densities of kT ~ 2.2 keV and n e ~ 
280 cm -3 , respectively. The lastest Chandra data indicate (Park et al. 2006) that 
the "fast shock" and the "slow shock" of the model are becoming less distinguishable, 
as the overall shock front is now entering the main body of the inner ring, that is, 
the electron temperature of the soft component is increasing and that of the hard 
component is decreasing. Albeit rather speculative, the overall temperature might 
thus be "merging" onto an "average" temperature of kT ~ 1.5 keV ~ 1.8 x 10 7 K, 
Park et al. (2006), and intermediate electron densities. The range of densities 
expected for this "average" shock is in very good agreement with that implied from 
the IR observations. 



4.2.2. The Infrared-to-X-ray Flux Ratio (IRX) 

An important diagnostic of a dusty plasma is the infrared-to-X-ray (IRX) flux 
ratio (Dwek et al. 1987b). For a given dust-to-gas mass ratio and grain size dis- 
tribution, the IRX ratio is defined as the infrared cooling to X-ray cooling in the 
0.2-4.0 keV band, and is given by: 

, s n. n d AJT ) 
n e n H A x (T e ) 

fim H Z d Arf(T) 
(m d ) A x (T e ) 
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were Z d = n d (m d ) / \xnuTnu is the dust-to-gas mass ratio, /i is the mean atomic 
weight of the gas, and (m d ) the mass of a dust particle, averaged over the grain size 
distribution. For a given value of Z d , the IRX ratio is only a function of plasma 
temperature, and ranges from a value of about 10 for T e = 10 6 K to a value of ~ 400 
for T e = 10 s K. For young supernova remnants, Dwek et al. (1987b) found that 
the IRX ratio is significantly larger than unity for 7 of the 9 remnants considered in 
their paper, the other two having only an upper limit on their IR emission. 

The observed IRX ratio can be obtained from the X-ray and IR fluxes from the 
SN. From the January 2005 Chandra data, we estimate the X-ray flux in the 0.2- 
4.0 keV X-ray band to be 7.24 x 10~ 12 erg cm~ 2 s _1 after correcting for interstellar 
absorption by an H-column density of Njj = 2.35 x 10 21 cm~ 2 . 

Note that these values are based on the two-temperature model as used in Park 
et al. (2004). Fractional contributions from the soft (kT e =0.3 keV) component in the 
total flux is ~70% for the unabsorbed flux (and thus for Lx)- Thus, a contribution 
from each component in the 0.2-4 keV X-ray flux seems to be significant rather than 
being dominated by one of them. 

The total IR flux is Fir = 7.7 x 10~ 12 erg cm -2 s _1 which leads to an IRX 
ratio ~ 1. This value is lower than the values reported in Paper I, in which IRX = 
6 for the decelerated slow shock component, and IRX = 3 for the blast wave shock 
in the two-temperature model. It is much lower than the theoretical value of ~ 10 2 , 
expected for a T e rs 1.8 x 10 7 K plasma, which is observed in the young remnants 
Tycho and Cas A. Other (mostly older) SNRs with measurable IR emission show 
somewhat lower IRX ratios, but not as low as SNR 1987A. 

Several effects could be the cause for this very low value of the IRX ratio in 
SNR 1987A. First, in remnants the IRX ratio was calculated by Dwek et al. (1987a) 
for an interstellar dust-to-gas mass ratio of 0.0075. The SNR 1987A blast wave is 
expanding into the circumstellar shell of its progenitor star, which, a priori, is not 
expected to have an interstellar dust-to-gas mass ratio. Moreover, when estimating 
the depletion of elements onto dust in the ER, we should compare the expected 
dust abundances in the ER with the maximum available abundances for the LMC. 
General LMC abundances exist for B stars (Rolleston et al. 1996) and ISM (Welty 
et al. 1999) and, although controversy still remains (Korn et al. 2002), the LMC 
metallicity is usually assumed to be 0.5 — 0.7 solar. Assuming that the fraction of 
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metals locked up in LMC dust is the same as in the local ISM, and that the ER has 
the same metallicity as the LMC, then the IRX ratio in the ER should be about 
0.5 — 0.7 times that expected from Supernova Remnants in the Milky Way, still 
significantly larger than implied from the observations. The extremely low value of 
the IRX ratio may therefore be due to a deficiency in the abundance of the dust, 
compared to interstellar values, which may reflect the low condensation efficiency of 
the dust in the circumstellar envelope. 

Second, the dust deficiency could be the result of grain destruction by thermal 
sputtering in the hot gas. The sputtering lifetime, T spu t, in a plasma with tempera- 
tures above ~ 10 6 K is about (Dwek et al. 1996): 

r sput » 3 x 10 5 a ^ m ] yr (10) 
n(cm d ) 

where n is the density of nucleons in the gas. The X-ray emitting gas is highly 
ionized, and we will assume that its density is that required to heat the dust to the 
observed range of temperatures, that is, n w 300 to 1400 cm" 3 . Grain destruction is 
important when the sputtering lifetime is about equal to the age of the shocked gas, 
which we take to be ~1 yr. The low IRX ratio can therefore be attributed to the 
effects of grain destruction if the dust particles had initial radii between 10 and 50 
A. So attributing the small IRX ratio to the effect of grain destruction in the hot 
plasma requires that only small grains had formed in the presupernova phase of the 
evolution of the progenitor star. 

The low IRX ratio shows that IR emission from collisionally heated dust is not 
the dominant coolant of the shocked gas. Its lower than expected value suggests a 
dust-to-gas mass ratio in the ER that is only a few percent of its interstellar value. 
The puzzle of the low dust abundance is greater if, in fact, the IR emission arises 
from dust that is not embedded in the X-ray emitting gas, but from dust that resides 
in the UV/optical knots instead. 



4.3. The Dust Heating mechanism 



The possibility that the initial size of the dust grains swept up by the shock is 
small suggests that the temperature of the grains may not be at the equilibrium value 



-18- 



but may fluctuate due to the stochastic nature of the heating and cooling. Very small 
grains will be stochastically heated if the energy deposited in the grain in a single 
collision is large compared to its internal energy content, and if its cooling time via 
IR emission is shorter than the time between subsequent collisions. Assuming that 
the electron and ion temperature is instantaneously equilibrated behind the shock, 
the mean thermal energy of the electrons will be 2.6 keV at a postshock temperature 
of 2 x 10 7 K. Using the electron ranges given by Iskef, Cunningham, & Watt (1983) 
we get that the energy, AE, deposited in a dust grain of radius a(/im) and density 
p(g cm -3 ) by electrons with energies E > 370 eV is given by: 

AE(erg) = 5.5 x 1(T 7 ^ , 492 (11) 

which for p = 3 g cm" 3 , a = 0.0050 jum, and E = 2.6 keV gives an energy deposition 
AE = 1.7 x 10~ 10 erg. This value is small compared to the internal energy of the 
grain, which is about 6 x 10~ 10 erg at the equilibrium grain temperature of 170 K. 
The average thermal speed of the electrons is about 3 x 10 9 cm s^ 1 , giving an average 
time between electronic collisions of about 1 sec for an electron density of 300 cm -3 
and a grain radius of 0.0050 /iin. This is shorter than the grain's cooling time at 
170 K, which is about 10 s (Dwek 1986), suggesting that the dust can maintain an 
equilibrium temperature at that value. 

Collisions with protons can deposit a significantly larger amount of energy in 
the dust. Typical proton energies in the shocked gas are ~ 5 keV, and the stopping 
power for protons with that energy is about 240 MeV cm 2 g _1 (NIST tabulated 
values: http://physics.nist.gov/cgi-bin/Star). The energy deposited in a 0.0050 /im 
radius dust particle is then about 6 x 10 -10 erg, which is about equal to the internal 
energy of the dust grain at 170 K. So collisions with protons and heavier nuclei can be 
neglected if they are in thermal equilibrium with the electrons, because of their lower 
collision rate. However, if the electron and ion temperatures are not equilibrated 
behind the shock, then the dust heating rate will be dominated by collisions with 
the protons. The time between successive collisions will be about 40 s for a proton 
density of 300 cm -3 . The grains should therefore cool to a temperature of about 
100 K before another collision will take place. However, the data does not support 
such a broad range of dust temperatures, suggesting that the gas density could be 
somewhat higher (by a factor of ~ 3), in which case the dust will cool only to a 
temperature of ~ 140 K, which may be more consistent with the observations. 
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All these scenarios support the idea that the dust temperature in the X-ray emit- 
ting gas does not fluctuate wildly about the equilibrium dust temperature, which can 
still provides strong constraints on the density and the equilibration of the electron 
and ion temperatures in the postshock gas. 

4.4. Dust in the Dense Knots of the Equatorial Ring 

The UVO light emitting knots discovered with the HST resemble a string of 
beads uniformly distributed along the ER. Figure 14 and Figure 15 depict the map 
of the dust optical depth overlayed with contour levels of the HST emission obtained 
on Dec 15, 2004 (day 6502), the closest to our 11.7 /im observations. The data look 
very similar, but the IR emission seems to emanate from a somewhat wider region 
than the optical emission, an effect that cannot be entirely accounted for by the 
lower resolution of the IR data. 

Nevertheless the good correlation between the IR emission maps and the HST 
image, suggest that a significant fraction, if not most, of the mid-IR emission may be 
emanating from the knots. The physical conditions of a particular knot (Spot 1 on 
the ER) have been modeled in detail by Pun et al. (2002), from the analysis of the 
UV/optical line emission detected by the HST Space telescope Imaging Spectrograph 
(STIS). They found that the UV fluxes could be fit with a model consisting of two 
shocks with velocities of v s = 135 and 250 km s -1 expanding into preshock densities 
of n = 3.3xl0 4 and 10 4 cm -3 , respectively. The postshock temperatures behind 
the slow and fast shocks are T s = 4.5 x 10 5 , and 1.5xl0 6 K, with gas cooling rates 
of A(T) ~ 10~ 22 erg cm 3 s _1 . The cooling time of the shocked gas is given by 
tcooi = kT s /nA(T s ), and is equal to w 0.1 and 2 yr for the slow and fast shock, 
respectively (Pun et al. 2002). The thickness of the shock front is therefore ~ 4x 10 13 
and 8xl0 14 cm for the slow and fast shocks, respectively, both significantly smaller 
than a typical radius of the knot, which is about 2 x 10 16 cm. So most of the dust 
in the knot resides in the unshocked gas and is heated by the radiation emitted from 
the cooling shocked gas. 
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The radiative energy density seen by the dust is approximately given by: 

c 

~ fc n okT s — 
c 

where £ coo ; = t coo \ v s , and / c is the compression factor of the gas in the postshock 
region. For T s = 10 6 K, no = 10 4 cm -3 , and a shock velocity of 200 km s" 1 we get 
that 

U rad « 2 x 10~ 9 x / c erg cm" 3 (13) 

Pun et al. (2002) find that compression factors can be as large as ~ 550, giving 
radiation densities of ~ 10~ 6 erg cm -3 throughout the knot. The energy density 
of the local interstellar radiation field is about 3 x 10~ 12 erg cm" 3 . Silicate dust 
particles immersed in this field achieve equilibrium dust temperatures of about 15 K 
(Zubko, Dwek, & Arendt 2004). The energy density in the knot is therefore higher 
by a factor of ~ 3 x 10 5 than that of the local interstellar radiation field, and the 
average dust temperature should therefore be higher by a factor of ~ 8.3 for a A" 2 

1 /6 

dust emissivity law (T d oc U rad ). This gives a typical dust temperature of ~ 125 K, 
in reasonable agreement with the observed average. 

The total mass of radiating dust was found to be ~ 10~ 6 M . The typical mass 
of gas in a knot of radius r = 2 x 10 16 cm, and density no = 10 4 cm" 3 is ~ 10~ 4 M & . 
If M ~ 20 is the number of knots in the ER, then the dust-to-gas mass ratio is 
~ 10~ 6 / (A/" x 10~ 4 ) ~ 5 x 10~ 4 , or approximately a factor of 10 less than the average 
dust-to-gas mass ratio in the local interstellar medium. 

The low abundance of dust in the knots could be explained if the dust is ef- 
ficiently destroyed in the shocked gas and if the transmitted shocks have already 
traversed most of the volume of the knots. Calculations presented by Jones (2004) 
show that about 49% of the silicate grains swept up by 200 km s" 1 shocks expanding 
into a medium with a preshocked density of 0.25 cm" 3 are destroyed. This fraction 
could be significantly higher for the densities encountered by the shocks traversing 
the knots. This scenario predicts that the IR emission from the knots was higher 
in the past, contrary to observed IR light curves (see below). Therefore, if the IR 
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emission emanates from the knots, the low dust abundance must reflect the initial 
dust abundance in these objects. 

5. THE LIGHT CURVES 

The light curves at 10 and 20 fim are shown in Figure 16: the absolute flux cali- 
bration have been made using Cohen et al. (1992) 0-magnitude fluxes [Fo(10.0 //m) 
= 35.24 mJy, F„(11.7 /im) = 28.57 mJy, and F (18.30 //m) = 10.25 mJy]. It can be 
seen from Figure 16 that the flux arising from the ejecta at 10 /im declines expo- 
nentially from day 2200 through day 4200 (the ISOCAM and OSCIR observations) 
until day 6000, at a rate of ~ 0.32 mag y~ x . This would imply that the observa- 
tions on day 4200 are not dominated by the ring emission. It is likely that the ring 
emission started around day 4000, at roughly the epoch when the first optical spot 
was discovered (Pun et al. 1997), and in good agreement with Fig. 5 of Park et al. 
(2002) which presents ATCA and Chandra/RO SAT data. The significant increase 
of the fluxes reported in the present paper compared to the previous fluxes (Fig. 16) 
is consistent with the soft X-ray flux increase observed in the last set of data (Park 
et al. 2005b). This clearly shows that "something" has happened at around day 
~ 6000. 

This is also manifest in the radio light curve at 843Mhz from MOST. The earlier 
observations by (Ball et al. 2001)from 3000 to 6000 days show a steadily increasing 
flux together with a steadily increasing rate of change. Near day 6000 however 
(Hunstead 2006) the rate of increase undergoes a more abrupt change upwards, and 
is therefore clearly associated with the changes seen at other wavelengths. Figure 16 
also shows the Spitzer data point resulting from integrating the IRS data within the 
appropriate bands. Including these data points shows that the trend of the increase 
of the flux at 11.7 fim since after day ~6000 is better fit by a "linear" function, while 
the increase in magnitude vs. time is not linear. Although we do not make strong 
conclusions, this could indicate that the shock is travelling through structures with 
cross sections (as seen from the center of the explosion) that do not increase with 
radial distance, like short cylindrical clouds. Sputtering time scales for all but the 
smallest grains are too long (10-100 years) to affect the IR emission, and the plasma 
cooling time due to the IR emission is also too long (~ 100 years) to affect the time 
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variation in the IR emission. 

The 13 cm radio emission light curve (http:/ /www.atnf.csiro.au/research/SN1987A/) 
is also shown in Figure 16. No data are currently available after day 6244 so it is 
premature to discuss the associated 13 cm evolution. It would be surprising if the 
rate of increase does not change upwards. While the sudden increase in the X-ray 
light curve at day around 3700 has been interpreted as the encounter of the shock 
front with the first protrusions of the ER (Park et al. 2004), the "jump" after day 
around 6000 could be the sign of the shock reaching the main body of the ER (Park 
et al. 2005b). These authors also show that the light curve of the hard X-rays (3-10 
keV) is much flatter than one corresponding to the soft X-rays, and similar to the 
radio light curves, and they argue that it is likely that the hard X-ray emission comes 
from the fast reverse shock rather than the decelerating forward shock, just like the 
radio emission (Manchester et al. 2005). We note that the reverse shock origin for 
the radio emission is one possible interpretation for explaining the inconsistency of 
the radio images with the IR and the soft X-rays images. It is not our intention to 
discuss evolution of radio fluxes and their relationship to other changes reported in 
this paper, but we note that the log-log plot in Figure 16 hides a significant change 
in the rate of increase of the 13 cm at 5000 days. The epoch of increase in the rate 
appears to be frequency dependent. 



5.1. The Ejecta 

An asymmetry in the profiles of optical emission lines that appeared at day 530 
showed that dust had condensed in the metal-rich ejecta of the supernova (Lucy et al. 
1991). Although it was discovered via spectroscopy, the presence of the dust could 
be easily inferred from the spectral energy distribution: as the dust thermalized the 
energy output, after day 1000 SN 1987A radiated mainly in the mid-infrared (Bouchet 
et al. 1991). Although the presence of the dust emission from the condensates in the 
ejecta was reported in Paper I, we do not detect them in the present observations. 
The possibility that this was due to the fact that the present observations were 
achieved with the narrow Si5 filter instead of the broad N filter was investigated, 
for the occurence of the [Ne II] 12.8 /iin line. Indeed, this line corresponds to a fine 
structure transition in the ground state of Ne II so the temperature to excite the 
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upper level does not need to be high if there is Ne II. However, Paper I reports a 
dust temperature of 90 K < T < 100 K in the ejecta, and there is little evidence 
for X-ray emission from the ejecta. Furthermore, the X-ray emitting gas would be 
ionized to a much greater extent than to produce Nell. It is thus most likely that 
Ne II is coming from warmer regions near the X-rays which are also responsible for 
the other HST and Spitzer lines. Why then is the ejecta not detected in our last 
observations? 

In Paper I, the N band background sigma was 0.033 mJy/pix. We used a 12 
pixels area to integrate the central source, and then the 3-cr detection limit was 0.34 
mJy which is about the 0.32 mJy measurement reported for the central source in 
Paper I. In the present data, although the source is observed with a better signal- 
to-noise ratio, the background at 11.7 fim is affected by a higher noise, with a 3-<r 
background value of 0.47 mJy. It transpires that we could not have detected a source 
at the flux level expected from the radioactive decay based light curve (eg. ~ 0.3 
mJy). It appears then that the only detection of the ejecta at this late stage is 
achieved in the H band (see Figure 7). The origin of this emission is most likely due 
to line emission, although we note that the limiting flux for any continuum emitter 
at the center of SN 1987A, in the wavelength range 2900-9650 A reported recently 
by Graves et al. (2005) (4.3 mJy in the I filter) is much above our detection at H 
(0.11 mJy). Therefore, our near IR results are compatible with HST observations 
and do not exclude contribution from a continuum. 

Supernovae are known contributors to interstellar dust. The presence of iso- 
topic anomalies in meterorites (Clayton & Nittler 2004) and observations of Cas A 
(Douvion et al. 1999; Arendt, Dwek, & Moseley 1999) and SNR 1987A (Paper I 
and references therein) provide direct evidence for the formation of dust in SNe. 
However, the relative importance of SNe compared to quiescent outflows from AGB 
stars in the production of interstellar dust grains is still unclear (Jones 2004; Dwek 
1998; Tielens 1998). So far, the dust masses in the the ejecta have been determined 
observationally for two SNe, SN1987A (Lucy et al. 1989) and SN1999em (Elmhamdi 
et al. 2003). These masses were ~ 3.10~ 4 M Q and ~ 10 -4 M Q . In both cases 
the authors note that these values could be much higher if the dust exists in opaque 
clumps which may be the case. At face value however they are much less than the 
0.1 - 1.0 M & required to make SNe the dominant source of interstellar dust particles. 
In future one must find a means of quantifying the effects of clumping. Observations 
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of the handful of other supernovae in which an IR excess is interpreted as dust form- 
ing in the ejecta (SN 1979C and SN 1985L, and probably SN 1980K) do not allow 
an estimate of the mass of dust. As for the young Galactic SNR which have been 
observed by IRAS and ISO (Cas A, Kepler and Tycho) the dust mass deduced is 
only KT 7 - 10~ 3 M (Lagage et al. 1996; ?), also many orders of magnitude lower 
than the solar mass quantities predicted. 

Paper I discusses the role of supernovae in dust production. We showed there 
that the dust which condensed in the ejecta of SN 1987A has survived 16 years since 
outburst, and was still radiating the energy released by the radioactive decay of 44 Ti 
at the expected level. Unfortunately, we could not accurately estimate the mass of 
the dust, and the observations reported in this paper do not allow it either, in the 
absence of detection of the ejecta. Thus, if we consider SN 1987A as an archetype, our 
data can neither support nor rule out the hypothesis that supernovae are significant 
sources for dust production. 



6. CONCLUSIONS 

We have presented mid-IR images of SNR 1987A obtained with T-ReCS on 
the Gemini South telescope on day 6526 at 11.7 and 18.3 jum and with IRAC (5.8 
and 8 /im) on day 6130 and MIPS (24 /im) on day 6187 onboard Spitzer, together 
with 3-37 /im spectroscopic observations of the remnant obtained with IRS on day 
6190 at the same observatory. The imaging observations (Figure 1) show that the 
mid-IR emission arises from the dust in the equatorial ring (ER) heated up by the 
interaction of the SN blast wave with its circumstellar medium. Several theoretical 
models predicted the presence of dust in the CSE of SN 1987A which was produced 
in the winds of the supergiant phase. The location of the IR emission rules out 
the possibility that the dust condensed out in the SN ejecta, strongly suggesting a 
circumstellar origin instead. 

The 3-37 /im spectrum (Figures 10) shows that the emission arises from a 
population of astronomical silicate particles. Temperature maps show that the dust 
temperature is fairly uniform in the ER and about 166l*2 K, with total dust masses 
of ~ 2.6±H x 10~ 6 M Q . 
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A comparison with Chandra and HST observations show an equally good correla- 
tion between the 11.7 /im IR and the X-ray and UV-optical images of the supernova. 
Because of the limited angular resolution in the IR we cannot determine the location 
or heating mechanism of the radiating dust. 

The dust could be residing in the hot ~ 10 7 K gas and collisionally heated by the 
X-ray emitting plasma. The dust temperature is then an excellent diagnostic of the 
electron density, giving a value of ~ 300-1400 cm -3 , similar to the value suggested 
by the Chandra observations. 

Comparison of the IR and X-ray fluxes suggests that the dust is depleted by a 
factor of ~ 30 in the X-ray emitting gas, compared to its value in the local interstellar 
medium of the LMC. This low value could be due to its destruction by thermal 
sputtering in the shocked gas, requiring the initial grain radii to be below ~ 50 A. 

Alternatively, the dust could be residing in the UV-optical emitting knots in 
the ER and radiatively heated by the cooling gas that was excited by the shocks 
propagating through these knots. A simple calculation for a particular knot shows 
that the radiative energy density can heat the dust to typical temperatures of about 
125 K, similar to those inferred from the IR observations. A comparison with the 
mass of the knots shows that the dust-to-gas mass ratio in the knots is lower by a 
factor of ~ 10 compared to its value in the ISM of the LMC. This low abundance 
reflects the low condensation efficiency of the dust in the outflow of the progenitor 
star. 

We stress that in order to assess the role of SNe in the production of dust in 
the Universe, it is clearly important to measure the presence of dust that survives 
into the formation of the remnant, and for this, mid-IR and sub- mm observations 
are critical. 
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Table 1. Observed Fluxes From SN1987A 1 



Day 2 Instrument 5.8 /im 8.0 /im 10.4 /zm 11.7 /iin 18.3 /im 24 /iin 

6067 T-ReCS ••• ••• 9.9±1.5 

6125 T-ReCS • • • • • • • • • • • • < 50.6 

6130 IRAC 1.79±0.06 4.98±0.16 

6190 IRS 3 ••• ••• ••• 13.3±0.3 31.8±1.0 

6184 MIPS ••• ••• ••• ••• ••• 29.8± 1.7 

6487 IRAC 2.44±0.19 7.48±0.19 

6526 T-ReCS ••• ••• ••• 18.4±1.2 53.4±9 



fluxes in units of mJy. 

2 Time elapsed since the explosion. 

3 Spitzer IRS data integrated between the filter half-power wavelengths. 



Table 2. Measured IR Lines from SN1987A at Day 6190 



Line Center (//m) 


Element 


FWHM (km s- 1 ) 


Line Flux (W cm 2 ) 


S/N 


N 


12.831 ±0.001 


[Ne II] 


442 


2.392 ±0.18510" 


-21 


20 


1 


14.337 ±0.004 


[Ne V] 


960 


7.153 ± 2.04310" 


-22 


4 


2 


15.577 ±0.001 


[Ne III] 


336 


1.038 ±0.14110" 


-21 


13 


1 


24.340 ± 0.002 


[Ne V] 


374 


4.819 ± 1.26510" 


-22 


9 


3 


25.833 ± 0.004 


[F IV] 


472 


4.613 ± 1. 16310 


-22 


7 


4 


25.904 ±0.001 


[0 IV] 


635 


2.016 ± 0.12210" 


-21 


24 


5 


26.005 ±0.002 


[Fe II] 


844 


1.683 ±0.14510" 


-21 


15 


5 


34.850 ± 0.008 


[Si II] 


603 


1.789 ±0.61610" 


-21 


4 


1 



^ine is strong in background. 

2 Line was only recognized by searching for counterpart to the 24.3 [im line. 
3 Line is faintly visible in 2-D spectral data. 

4 Line is not well-resolved. It could also be fit as the wing of a very broad line 
underneath the [O IV] and [Fe II] lines. 

5 Line is weak or absent in background. 
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Table 3. Global parameters from Silicate Dust Modeling 1 



Instrument 


Day 2 


Tdust (K) 


Lir (erg s 1) 


M dust (M ) 


T-ReCS 
Spitzer 
T-ReCS 


6070 
6190 
6530 


180+ 2 ° 
180 +2 ° 
1661J1 


9 ± 3 x 10 35 
1.6 ±0.3 x 10 36 
2.3 ±0.4 x 10 36 


1 - 4 x 10~ 6 
0.7-1.7 x 10~ 6 
1.2-4.6 x 10~ 6 



1 The dust mass for day 6070 has been recalculated from Paper I with the value of 
k adapted to the silicate case. 

2 Time elapsed since the explosion 
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Figure 1 
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Fig. 1. — (a) SN 1987A seen with T-ReCS at day 6526 in the Si5 narrow band filter 
(11.7 /im) and (b) at day 6552 in the Qa filter (18.3 /xm); (c) temperature map 
assuming pure black body emission and Mathis (1990) extinction law; (d) Opacity 
map resulting from the same algorithm. These images are smoothed 2 pixels (0.18 
arcsec) . 
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Figure 2 
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Fig. 2. — Comparison of the N filter (10 //m) emission at day 6067 with the Si5 filter 
(11.7 /im) emission at day 6526: (a) flux calibrated difference of the two images; (b) 
that difference divided by the 2003 image; contours are drawn at [-50,150,250,350,700] 
percent changes; (c) image with contours from the 2003 image; (d) image with contours 
from the 2005 image. All the images are smoothed 2 pixels (0.18 arcsec). 
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Fig. 3. — Overlays of the contour images obtained with ACIS on Jan. 9-13, 2005 
at the Chandra X-ray Observatory superimposed on the T-ReCS Si5 (a) and Qa (b) 
image. The T-ReCS image has been smoothed 2 pixels (0.18 arcsec). 
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Figure 4 




Fig. 4. — Upper: contours from the ATCA image obtained on Jul. 31, 2003 in the 
12 mm band (16-26 GHz) superimposed on the T-ReCS Si5 (a) and Qa (b) images; 
lower: ATCA image at same frequencies obtained on May 5, 2004 (ATNF Web page) 
with contours from T-ReCS at 11.7 /im (c) and 18.3 /xm (d). The T-ReCS image 
have been smoothed 2 pixels (0.18 arcsec). 
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Figure 5 
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Fig. 5. — The Spitzer IRS spectrum. The blue and green symbols are from the short 
low resolution mode; the yellow and red are median binned data from the short and 
long high resolution modes. The black dots are aperture photometry data points 
from IRAC and MIPS. 
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Figure 6 




Fig. 6. — Lines emission from SN 1987A observed by the IRS on theSpitzer Space 
Telescope (see text). 
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Figure 7 
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Fig. 7. — Top: SN 1987A seen with ISPI in the J, H, and K band at the 4-m Blanco 
telescope at CTIO; the images have been deconvolved using the Multiscale Maximum 
Entropy Method from Pantin & Starck (1996). Bottom: near-IR images obtained 
with HST on Nov. 11, 2005, in the F110W (0.8 - 1.4 H. F160W (1.4 - 1.8 /mi), 
and F205W (1.75 - 2.35 /im) filters; the T-ReCS 11.7 fim image is also shown at the 
same scale. 



-40- 



Figure 8 
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Fig. 8. — The values of the dust mass absorption coefficient, for the three types 
of grains considered (references are given in the text). For small dust particles with 
radii < A, the value of k(X) is independent of grain radius . 
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Figure 9 




A(/-im) 

Fig. 9. — The T-ReCS fluxes at 11.7 and 18.3 /im fitted by our modelling of the dust 
continuum for 3 different dust compositions: silicate (blue line), carbon (green line), 
and graphite (red line) and in the simple black body case (black line) . Note that the 
T-ReCS data alone cannot determine the dust composition. 
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Figure 10 
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Fig. 10. — The Spitzer IRS spectrum fitted with a silicate dust model. Green cir- 
cles: Spitzer flux measurements. This model has been scaled to fit the T-ReCS flux 
measurements (red squares). Note that the mid-IR emission has been brightening 
significantly between both sets of observations. The Spitzer IRS data unambiguously 
show that silicate are the major dust component. The derived parameters of the fit 
are: T = 180±|g K, M dust = 0.7 - 1.7 x 1O~ 6 M , and L IR = 1.6 ± 0.3 x 10 36 erg s -1 . 




Fig. 11. — Residuals from Spitzer observations fitted with our silicate dust mo 
The green curve shows emission from amorphous carbon dust at a T = 750 K. 
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Figure 12 




Fig. 12. — Temperature (a) and optical depth r(11.7 /im) (b) maps derived for the 
T-ReCS observations for a silicate dust composition. 




Fig. 13. — Contours of silicate dust temperature collisionally heated by a hot gas as 
a function of plasma temperature and electron density. 
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Figure 14 




Fig. 14. — (a) Montage of the mid-IR 11.7 /im image of the ER (red) and the HST 
F625W image obtained on day 6502 (yellow); (b) same as (a) with the HST image 
(yellow) convolved to match the resolution of the T-ReCS 11.7 fim image (red); (c) 
montage of the Qa image from T-ReCS (red) with the HST F625W image (yellow); 
(d) same as (c) with the HST image convolved to match the T-ReCS Qa angular 
resolution. 
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Figure 15 




Fig. 15. — (a): Image of the ER obtained through the 11.7 /im filter with T-ReCS 
with contours of the HST UVO image obtained on day 6502 superimposed; (b): 
same with the 11.7 /im image deconvolved with the "Multiscale Maximum Entropy 
Method" from Pantin & Starck (1996). 
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Figure 16 
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Fig. 16. — Light curves of SN 1987A for the indicated frequencies; the X-ray 
data points before day ~ 4000 are from the ROSAT data (Hasinger,Aschenbach, 
B., & Truemper 1996); ATCA and ROSAT/ACIS fluxes are plotted in ar- 
bitrary logaritmic scales; ATCA data points are from the ATNF web page 
(http:/ /www. atnf.csiro.au/research/SN1987A/) and ACIS data are from Park et al. 
(2005b). The upper point in the Spitzer/IRS data is integrated in the 18.3 /im band- 
pass, and the lower point is integrated in the 11.7 //m bandpass. The point called 
'Ejecta' derives from the point source near the center of the image reported in Paper 
I. After day 4000 all data points refer to the inner ER. 
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Fig. 17. — Mid-IR images of SN 1987A from Spitzer, compared to a near-IR im- 
age from the Hubble Space Telescope. Left to right, top to bottom: (a) Hubble 
ACS/HRC F814W (range = [30, 3.2e4] electrons, pixel = 0.05"); (b) Spitzer IRAC 
3.6 fan (range=[-0.2, 4] MJy/sr, pixel = 1.2"); (c) Spitzer IRAC 4.5 /im (range=[0.2, 
5] MJy/sr, pixel = 1.2"); (d) Spitzer IRAC 5.8 /im (range=[1.8, 8] MJy/sr, pixel = 
1.2"); (e) Spitzer IRAC 8 /im (range=[2, 20] MJy/sr, pixel = 1.2"); (f) Spitzer MIPS 
24 fim (range=[18, 45] MJy/sr, pixel = 2.45"). The bright stars in the Hubble image 
NW and SE of SN 1987A (Star 2 and Star 3) appear to be brighter than the SN 
emission at 3.6 /im, but rapidly fade as the SN brightens at longer wavelengths. The 
large ring in the MIPS image is the first diffraction ring of the PSF, which is visible 
for any bright unresolved source. All pictures are oriented with N upward and E 
to the left. All span a 35.2" x 34.6" field. The Hubble image is a preview image 
from the Multimission Archive at Space Telecope (MAST). The Spitzer images are 
post-BCD images from the science data archive at the Spitzer Science Center. 
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Fig. 18. — IRS slit positions overlaid on the Spitzer IRAC 8 /im image (range = [0,15] 
MJy/sr). This image spans a 82.8' square field. The slit positions indicated by the 
long narrow red/yellow rectangles are the IRS short-lo (SL) lst/2nd order (7.4-14.2 
/ 5.3-8.5 jum) pointings. The short narrow blue rectangles indicate the short-hi (SH) 
pointings. The short wide green rectangles indicate the long-hi (LH) pointings. For 
each IRS module the target is positioned at locations 1/3 and 2/3 of the distance 
along the slit length. 



